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Abstract. We calculate ion distributions near a planar oil- water interface within 
non-linear Poisson-Boltzmann theory, taking into account the Born self-energy of the 
ions in the two media. For unequal self-energies of cations and anions, a spontaneous 
charge separation is found such that the water and oil phase become oppositely charged, 
in slabs with a typical thickness of the Debye screening length in the two media. 
From the analytical solutions, the corresponding interfacial charge density and the 
contribution to the interfacial tension is derived, together with an estimate for the 
Yukawa-potential between two spherical water droplets in oil. The parameter regime 
is explored where the plasma coupling parameter exceeds the crystallization threshold, 
i.e. where the droplets are expected to form crystalline structures due to a strong 
Yukawa repulsion, as recently observed experimentally. Extensions of the theory that 
we discuss briefly include numerical calculations on spherical water droplets in oil, and 
analytical calculations of the linear PB-equation for a finite oil-water interfacial width. 
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1. Introduction 

The making and breaking of oil-water emulsions is not only a problem of extreme 
importance in chemical, oil, pharmaceutical, food, and cosmetics industries, but is also 
a scientifically fascinating topic. It is well-known that the intrinsic tendency of oil and 
water to demix can be slowed-down or delayed by adding surfactants or colloidal particles 
to the mixture. These additives strongly adsorb to the oil-water interface, which lowers 
the interfacial tension (the main driving force for demixing) and/or provides a kinetic 
barrier that prevents droplet coalescence [HE]. Recently however, it was observed that 
emulsions of water droplets dispersed in somewhat polar oils can be stable for a long 
time (so far for more than 18 months), without any additives [3]. Moreover, the water 
droplets, which are of micrometer dimensions in these experiments, spontaneously form 
crystalline structures with lattice spacings of the order of 5 — 15 /im. These observations 
suggest long-ranged electrostatic droplet-droplet repulsions due to a net water droplet 
charge stemming from a preferential uptake of ions from the oil [31 HI [5] . Theoretical 
calculations based on Poisson-Boltzmann theory for monovalent ions in the geometry of 
a planar water-oil interface showed that the order of magnitude of the charge separation 
process, caused by the different cationic and anionic Born self-energies in oil and water, 
is indeed sufficiently strong to explain the observed stability and crystallization [HE], 
at least qualitatively. These predictions are based on the assumption of a pairwise 
screened- Coulomb potential of the charged water droplets through the oil, and an 
explicit empirical crystallization condition based on simulations j6j [7] . 

In this paper we will explore the ionic charge separation at the planar oil-water 
interface and the resulting crystallization regime of water droplets dispersed in oil in 
the high- dimensional parameter space of salt concentration, dielectric constant of the 
oil, ionic sizes and self-energies, droplet size, and droplet concentration. The analytic 
solution admitted by the nonlinear Poisson-Boltzmann equation in the planar geometry 
[HI [9] allows for such a detailed exploration efficiently; the effects of droplet curvature 
were recently studied numerically [10]. We will find that crystallization of water droplets 
is only possible for sufficiently large droplet radii (> 100 nm), sufficiently large water 
content (volume fraction > 10 -3 ), and sufficiently (but not too) polar oils with dielectric 
constants between 4 and 10. In addition we will also give an explicit expression for the 
change of the interfacial tension between oil and water due to the ionic segregation over 
the two media. This change is found to be negative, proportional to the square-root of 
the ionic strength, and small on the scale of the bare oil-water tension, at least within 
the present model. We will discuss the sensitivity of these findings to details of the 
interfacial structure as described recently in [II] . 

2. Poisson-Boltzmann theory for a planar oil- water interface 

The liquids are considered as structureless homogeneous linear dielectric media, filling 
the two half spaces z < (water) and z > (oil), forming a flat interface at z — 0. The 
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relative dielectric constant is a step function, e(z) = e w (z < 0) and e(z) = e Q {z > 0). 
The grand potential functional per unit area of the variational density profiles p±{z) of 
the cations (+) and anions (— ) in units of k^T = 1/(3 is given by [51 fTUl [TT1 [12] 



with the self-consistent dimensionless electrostatic potential <f>(z, [p±]), and with the 
external potential acting on the ions 



representing the Born self-energy and an additional specific solvation energy g± (e.g. due 
to hydration, hydrogen bonding, local modification of dielectric constant [T3"l [T4"t IT5]). 
which is set either to or to 4 here in order to study specific effects. The salt 
concentration in bulk water is denoted by p w , e is the elementary charge, and a± is 
the ionic radius. Both <fi(z) and V±(z) are gauged to zero in the bulk water phase 
(z — > — oo). Figure [1] shows f± as a function of the oil dielectric constant e Q , and ionic 
radius a± (inset), for g± = 0. Typical values for the self-energies are 5 — 20 k^T, and 
self-energy differences are of the order of 1 — 10 k^T for ionic radii differing by, say, 1 
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Figure 1. The Born self-energy f± in units of k^T as a function of dielectric constant 
e , for several ionic radii a±, and as a function of the ionic radius at e Q = 7.5 (inset). 



Minimization of the functional (0Q) with respect to p±, together with the Poisson 
equation, yields 

p±(z) = Pw exp(-[3V±(z)T^z)), (3) 
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which can be rewritten as the Poisson-Boltzmann equation 

J? (A = I 4 sinh (0W) ,z < 0; ( . 

subject to the Neumann boundary conditions 
lim 4>'(z) = lim (j)'(z) = 0; 



-oo 



lime w 0'(z) = lime o 0'(z), (5) 

which are dictated by the vanishing of the electric displacement field in the bulk and 
its continuity at the interface. The Donnan potential of the bulk oil is <f>(z — > oo) = 
0d = §(/- — /+)) the inverse Debye screening lengths are given by k- x = ^8ire 2 p-Je^ksT, 
i G {o,w}, and the bulk salt concentration in oil satisfies p Q = p w exp(— -(/_ + /+))• 
The general solution of (j4]) is 

4 arctanh(C w e KwZ ) , z < 0; , . 

4 arctanh(C e- KoZ ) + D , z > 0, ^ ' 

with integration constants that follow from the boundary conditions (jSJ) as 

„ n + cosh 4p- — -D 
— f 

sinh 

1 + n cosh 4? - D 

C = 1 , (7) 

n sinh ^p- 



where n = e w K w /e K , and D = yn 2 + 2n cosh 4f- + 1. The derivation so far is 
equivalently presented in [8j [9]. A separation of charge is found near the interface 
for unequal self-energies, comprising a cloud of net charge in the water phase of typical 
width K" 1 and one in the oil phase of typical width k~ 1 , with a charge per area at the 

water side given in units of e by 

o 

= J dz(p + (z) - p_{z)) = ~ 8 / f W 1 _^2 ' ( 8 ) 

— oo 

the compensating charge resides at the oil side as the system is globally neutral. In 
section [3] we will assume that the surface charge density of a water droplet in oil is equal 
to cr w in order to calculate the pair potential between the droplets. 

3. Screened Coulomb potential between water droplets in oil 

We now consider iV oil-dispersed water droplets of radius a in a volume V, such that 
the typical droplet-droplet separation is R = (V/N)3 and the water volume fraction 
is x = ^7ra 3 N/V = ^ii(a/R) 3 . We assume that each water droplet has a charge 
Z = 47ta 2 o~ w , where a w follows from the nonlinear PB theory (jSJ), and reconsider a 
screened Coulomb (Yukawa) potential between two droplets 
rr . , {Zef ( e K ° a \ 2 e" v 
e Q \1 + K a/ r 
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with r the center-to-center separation between the droplets. According to simulation 
results, crystallization of such a Yukawa system occurs if V > 106 [61 [7J, with the 
coupling parameter Y defined by 

T = f3V Ynk (R)(l + k q R + ^!). (10) 

Although this condition has been confirmed only for point- Yukawa systems, it is 
expected to hold also for finite droplets as long as VYuk(2a) > 10 k^T and n Q a < 1, 
conditions that are easily met for micron-sized droplets, provided 4 < e D < 10 as can be 
seen in figure [3j 

The electrostatic contribution 7 to the interfacial tension can be calculated 
analytically by evaluating the functional ([1]) with the equilibrium profiles, minus the 
functional evaluated with the step profile p±(z < 0) = p w , p±(z > 0) = p Q [101 [TTJ [15]. 
After a tedious but straightforward calculation we find 

h = - 16pw 1 C % (c w - c a ) oc (ii) 

which is symmetric under exchanging w <-> o and where the proportionality to y / /v 
follows from and the fact that C w and C Q are independent of p w . Typical 

values of the excess interfacial tension are at most 0(pN /m), and are therefore 
vanishingly small compared to the value of the bare oil- water interfacial tension ((9(10 
mN/m)). 



4. Numerical results 



As standard parameters we use e Q = 7.5, p w = 10~ 3 M, which are close to the 
experimental values of [31 H]> g± = and (a + ,a_) = (0.36,0.3) nm. The surface 
charge density <r w of the interface appears to be strongly dependent on the oil dielectric 
constant, keeping the ionic radii fixed. Typically it is of the order of O(10 -4 ) elementary 
charges per nm 2 in the range 10 < e Q < 70, and decays rapidly to 10~ 7 nm -2 for e Q = 4, 
see figure [2](a). Numerical calculations in the spherical geometry [10] predict a smaller 
cr w for oil-in-water droplets, and a larger one for water-in-oil droplets in the regime 
4 < e Q < 10, as shown in figure [2]^a) for a radius a = 1 pm. The excess interfacial tension 
7 differs correspondingly for finite droplets, figure [2](b). The linearized PB theory that 
takes into account a finite interfacial width s ^ [H] between oil and water agrees 
quantitatively with the present results for <r w , indicating that the charge separation is 
hardly dependent on s = (9(10 -1 nm), even if the ions are effectively excluded in a 
band of several tenths of nanometers. On the other hand, above a certain crossover 
ionic strength, the same theory predicts a qualitatively different asymptotic behaviour 
of the excess interfacial tension, being positive and proportional to the ionic strength 
p w , figure EJ^b) . One can understand that as follows: whereas the charge separation is 
hardly affected by the interfacial width s, the adsorption of particles, and hence the 
interfacial tension, is sensitive to the effective exclusion in a region of width s. As the 
absolute value of the excess interfacial tension is vanishingly small compared to the bare 




Figure 2. (a) Interfacial charge density as a function of the oil dielectric constant, (b) 
Excess interfacial tension as a function of ionic strength, at e D = 7.5. If the interfacial 
width s — 0, the planar calculations clearly show a square root behaviour. The sign 
of 7 is negative. The dashed lines show the numerical results from the PB-theory in 
spherical geometry, for water droplets of radius a = 1 /xm (WO), and equal sized oil 
droplets in water (OW) [10]. The thin lines show the result of the planar linearized 
PB-theory with interfacial width s, for which 7 has a linear asymptote and positive 
sign in the high-p w limit [11] , if s =/= 0. 

oil-water interface, 7 ow = O(10 mN/m), 7 will be considered not to contribute to the 
stability of the syst 
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Figure 3. The Yukawa potential at contact and the Debye screening length in oil 
(inset) as a function of the oil dielectric constant e Q , at p w = 10~ 3 M, (a+,a_) = 
(0.36,0,3). The contact value VVuk(2a) > 1 k^T only for e > 4 and droplet radii 
a > 100 nm. On the other hand, the range of the Yukawa potential decays rapidly for 
e D > 10 (inset); for droplets with a radius a > 100 nm, the value of n a ^> 1 above 
e Q > 12. Therefore, significant long range Coulomb interactions are only expected 
between 4 < e Q < 10. 
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We expect that water droplets in oil can only form crystalline structures if long range 
interactions are present (k cl < 10), and if the repulsion is sufficiently strong. Figure 
E] shows the Yukawa potential ([9]) at contact, VV u k(2a), for several droplet sizes as a 
function of e Q , revealing that these two conditions already impose strong restrictions on 
the dielectric constant and size of the particles. The inset of figure [3] shows for instance 
that the screening length in oil decays from 10 — 100 fim at e G = 4 to 100 nm at e Q = 12, 
for p w = 1 mM, while the main part of figure [3] shows contact potentials that exceed 
10 k-gT for e Q — 6 — 12 provided the droplet size is in the micron regime. Thus one 
only expects crystallization to be possible for micron-sized water droplet in oils that are 
sufficiently polar (e Q > 4) to have enough charge (see also figure EJ^a)), but not too polar 
(e Q < 12) to have a long enough range k~ 1 of the repulsions. This regime of e D is in 
remarkably good agreement with the experimentally found regime 0,1!]. 

The actual parameter regime where droplets are expected to form crystalline 
structures, on the basis of T > 106, is depicted in figure HJ where the varied parameters 
are the ionic radii a + ,a_ (figure Hta)) and the volume fraction and droplet radius x, a 
(figure IHb)). The lines show the envelopes of the regimes T > 106 for all physically 
achievable ionic strengths p w . Variations of the external potential (5V±{z) by taking a 
non-zero g± are seen to result into a significantly modified crystallization regime. Other 
solvation effects, e.g. hydration, hydrogen bonding, local alignment of dipolar fluid 
molecules, are therefore expected to be important for a proper quantitative picture. We 
examined this for g± = 4 (figure HJ, which is considerable compared to \ f + — f_\ ~ 2 for 
the present parameters but small compared to f± ~ 12. Figure H(b) shows that micron- 
sized droplets, for our standard parameter set, tend to crystallize at x ~ 10~ 2 — 10 _1 , 
which is somewhat higher than the experimentally observed regime x ~ 10~ 3 — 10 -2 
[3|,|1]. A similarly relatively high theoretical value for the crystallization volume fraction 
was found in [10], where curvature effects were studied, with g± = 0. Combining these 
results with those of figure BJb) suggests that the experimental results may only be 
quantitatively explained by taking both curvature and specific solvation effects (g± ^ 0) 
into account. This problem is left for future studies. 

5. Concluding remarks 

By means of a non-linear Poisson-Boltzmann theory, that takes into account the Born 
self-energies of the ions [8j [9], and an specific solvation energy [T3J [TJJ [15], ionic 
equilibrium distributions near an oil-water interface were calculated. A difference in 
self-energy was found to lead to a Donnan potential and charge separation near the 
interface. The analytical expression for the charge of the medium per area of interface 
agrees roughly with numerical results of the Poisson-Boltzmann theory in spherical 
geometry [10] , although the planar charge density is somewhat underestimated compared 
to the value for micron-sized water droplets in oils with e Q < 10. The charge is in 
good quantitative agreement with a linearized Poisson-Boltzmann theory that takes 
into account a finite interfacial width s |11] , showing independence of the charge on 
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Figure 4. (a) The ionic radii (a+, a_) for which T > 106 for some physically achievable 
ionic strength p w . The droplets crystallize in the area between the curves and the 
horizontal axis (R/a = 10), and between the curves and the vertical axis (R/a — 5), 
respectively. The lines can be mirrored in the diagonal since there is a symmetry under 
a + <-> a_ , together with g + <-> g_ . The thin and dashed lines show the results if one 
of the ionic species has an additional solvation energy (independent of p w and e Q ) of 
g± = 4. The dielectric constant of the oil is e = 7.5 and droplet radius a = 1.5 /im. 
(b) The volume fractions of water x and droplet radii a for which T > 106 for some 
physically achievable ionic strength p w , for several combinations for g±. 



s. For slightly polar oils, i.e. 4 < e Q < 10, the typical charge of micron-sized water 
droplets ranges between 10 — 1000 elementary charges and the Debye screening lengths 
range between 0.1 < n a < 10, such that the particles repel each other electrostatically 
by many k^T even if their centers are several radii apart. On the basis of an empirical 
crystallization condition [7] we found a relatively narrow regime in the high dimensional 
parameter space, where water droplets in oil are expected to form crystalline structures, 
namely only when 4 < e Q < 10 (in striking agreement with experimental findings [3]), 
x > 0.001, a > 100 nm, and a sufficiently large difference between the ionic self-energies 
\f+ — f-\ > 1 k^T. The spontaneous charging of the water droplets might play an 
important role in the production and stabilization of emulsions with rather polar oils, 
since we have seen that by a judicious choice of the type of oil and salt ions, water 
droplets can be stabilized by these, sometimes surprisingly strong, electrostatic effects 
alone. 
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